Spin dynamics of heterometallic CryM wheels (M = Mn, Zn, Ni) probed by inelastic 

neutron scattering 



O 

o 

(N 
o 

Q 
o 



C/3 



o 

(N 
> 

O 



O 

o 



X 
S3 



R. Caciuffo and T. Guidi 

Dipartimento di Fisica ed Ingegneria dei Materiali e del Territorio, 
Universitd Politecnica delle Marche, Via Brecce Blanche, 1-60131 Ancona, Italy 

G. Amoretti, S. Carretta, E. Liviotti, and P. Santini 
Dtpartimento di Fisica, Universitd di Parma, F43100 Parma, Italy 

C. Mondelli 

Istituto Nazionale per la Fisica della Materia, and Institut Laue Langevin, 
Boite Postale 220 X, F-38042 Grenoble Cedex, France 

G. Timco, C. A. Muryn and R.E.P. Winpenny 

Department of Chemistry, The University of Manchester, 
Oxford Road, Manchester, MIS 9PL, United Kingdom 

Inelastic neutron scattering has been applied to the study of the spin dynamics of Cr-based 
antiferromagnetic octanuclear rings where a finite total spin of the ground state is obtained by 
substituting one Cr'^"'" ion (s — 3/2) with Zn (s — 0), Mn (s — 5/2) or Ni (s = 1) di-cations. Energy 
and intensity measurements for several intra-multiplet and inter-multiplet magnetic excitations allow 
us to determine the spin wavefunctions of the investigated clusters. Effects due to the mixing of 
different spin multiplets have been considered. Such effects proved to be important to correctly 
reproduce the energy and intensity of magnetic excitations in the neutron spectra. On the contrary 
to what is observed for the parent homonuclear Crg ring, the symmetry of the first excited spin 
states is such that anticrossing conditions with the ground state can be realized in the presence of an 
external magnetic field. Heterometallic CrrM wheels are therefore good candidates for macroscopic 
observations of quantum effects. 

PACS numbers: 75.50.Tt, 75.10.Jm, 75.40.Gb, 75.45.+j 



I. INTRODUCTION 

Magnetic wheels are polynuclear molecular clusters 
with a ring-shaped cyclic structure and a dominant an- 
tiferromagnetic (AF) coupling between nearest-neighbor 
ions. For an even number of spin centers, the ground 
state is a singlet and the excitation spectrum is charac- 
terized by rotational and spin- wave like bands Q, |^ 1^ . 
Heterometallic rings with can be obtained from 
an S = homonuclear ring by chemical substitution of 
one or two magnetic centers 0, |^ . Theoretical calcula- 
tions suggest that such systems could exhibit interesting 
quantum-coherence phenomena Q and are therefore 
of considerable interest. Indeed, the replacement of a 
magnetic ion in a cyclic structure allows one to modify 
the topology of the exchange interactions, that in turn 
plays a key role in determining the macroscopic behavior 
of the system HHIl^. 

Here we report the results of inelastic neutron scatter- 
ing (INS) experiments on heterometallic AF rings, and 
we show that the spin level sequence and dynamics are 
substantially modified with respect to those of the parent 
compound. The investigated wheels are derived from the 
spin-compensated neutral [Cr8F8(02CCMe3)i6] ring 
by substitution of one divalent cation (M = Mn, Zn, Ni) 
for a trivalent Cr ion Due to the difference between 
the Cr'^^ spin (s = 3/2) and the spin of the dication, the 



ground state of the so-obtained CryM wheels has a non- 
zero total spin (S = 1/2, 1 and 3/2 for M = Ni, Mn and 
Zn respectively). Recent investigations have shown that, 
within this CryM family, S = 3/2 rings exhibit magne- 
tocaloric effects that could be exploited below T ~ 2 K 
p^. whilst CrrNi could be a suitable candidate for the 
physical implementation of qubits 0. The suggested 
opportunities are linked to properties that are very sen- 
sitive to the spin energy spectrum and to the compo- 
sition of the spin wavefunctions. A reliable determina- 
tion of these quantities is therefore relevant, and INS is 
known to be the most appropriate technique to address 
this problem p_4 d5i. .l6nl7.J . 

Our INS experiments on CryM provide energy and in- 
tensity for several transitions between the anisotropy- 
split lowest spin multiplets. The analysis of the data 
leads to an accurate determination of the main features 
of the microscopic intra-cluster interactions, including 
nearest-neighbor isotropic exchange, zero-field splitting 
and dipole-dipole interaction parameters (inter-cluster 
interactions are very weak, of the order of 10 neV). The 
results obtained indicate that anticrossing conditions be- 
tween the ground state and excited low-lying energy lev- 
els with the same symmetry can be met in CryM, by 
applying a suitable magnetic field. This is at variance 
with the situation in the parent Crg ring, where ground 
state crossings involve states with different symmetry 



[l5j, as confirmed by the vanishingly small level repul- 
sion observed by heat capacity measurements The 
occurrence of anticrossing conditions suggests that CryM 
wheels are good candidates for the macroscopic observa- 
tion of quantum phenomena arising from S-mixing, such 
as the quantum fluctuation of the magnitude of the total 
spin d^J. 

A correct interpretation of the experimental spectra 
requires a generalization of the orientation-averaged INS 
cross-sections for polycrystals , as proposed by Wald- 
mann ,2Q |. 



II. EXPERIMENTAL DETAILS 

Deuterated micro-crystalline samples of (C2D5)2NH2 
[Cr7MF8(02(C5D9))i6 ] (M2+ = Ni, Mn, Zn) have been 
prepared according to a slightly modified literature pro- 
cedure 0i 0) by dissolving chromium fluoride tetrahy- 
drate (Aldrich) in a mixture of trimethyl-dg-acetic acid 
and diethyl- rfio-amine (98% deuterated, Aldrich) be- 
fore adding an excess of the second metal salt (nickel 
carbonate hydroxide tetrahydrate, manganese chloride 
tetrahydrate or basic zinc carbonate). The preparation 
of trimethyl-rfg-acetic acid startin g fro m acetone- rfg was 
adapted from standard methods [21|. All compounds 
were purified additionally on a silica gel column using 
toluene as the eluent Q and finally crystallized from a 
mixture of pentane/acetone by evaporation of the mix- 
ture of solvents at 313-318 K. X-ray diffraction analy- 
ses for all three compounds showed that crystals don't 
have any solvent molecules in the crystal lattice. The 
observed Bragg peaks could be indexed in the tetragonal 
PA space group, with lattice parameters at 100 K of a 
= b = 19.9598(2) A, c = 16.0609(2) A for CryNi; a =b 
= 19.9375(2) A, c = 16.1079(2) A for CryMn and a =b 
= 19.8973(2) A, c = 16.0227(2) A for CrrZn. The mag- 
netic rings have a diameter of about 10 A and the metal 
ions have the shape of an almost regular planar octagon 
(Figure P). 

INS measurements have been performed at the In- 
stitute Laue-Langevin in Grenoble (France), with the 
direct-geometry time-of- flight spectrometer IN5, using 
about 2 grams of each compound. A flat disk geom- 
etry has been chosen for the aluminum sample holder 
(1.1 mm thickness, 5 cm diameter), in order to reduce 
multiple-scattering background. Neutron spectra have 
been recorded with the sample kept at different temper- 
atures, between T = 2 K and T = 12 K, inside a standard 
liquid-'* He cryostat. Solid angle and detector efficiency 
calibrations have been performed using the spectrum of 
a vanadium metal sample. Neutron incident wavelengths 
A = 4 A , 5 A and 9 A were used, corresponding to in- 
strumental resolutions at the elastic peak of 0.147 meV, 
0.091 meV, and 0.019 meV respectively. The three sets 
of measurements allowed us to cover an energy-transfer 
range from hui ^ 0.04 meV to about 4 meV. The angular 
interval spanned by the detector banks corresponds to 




FIG. 1: The molecular structure of the CrrM (M = Mn, Ni, 
Zn) magnetic ring, with hydrogen atoms omitted for clarity. 
The A ions (A = 7/8 Cr and 1/8 M at each site) are repre- 
sented by large dark circles; fluorine positions are indicated by 
small dark circles, the oxygen and carbon atoms are large and 
small gray circles, respectively. The cation is found H-bonded 
inside the cavity. 



scattering vector amplitudes varying between Q ~ 0.55 
A~^ and 2.6 A~* at Hlu ^ I meV, and between Q ~ 0.85 
A-i and 2.2 A"* at Hlu - 3.7 meV (for A = 4 A ). 



III. SPIN HAMILTONIAN AND INS 
MAGNETIC CROSS SECTION 

The investigated systems are ensembles of indepen- 
dent identical magnetic units, each one described by the 
Hamiltonian operator: 

6 

if = ^ Jcr-CrSi ■ Si+i + Jcr-A/(S7 ' Sg + Sg • Si) + 
i=l 

8 8 
-I- ^ Si ■ D; ■Si+ ^ Si • Dy • Sj. (1) 

i=l i<j = l 

In the above equation we assume that sites i = 1 — 7 
of the octanuclear ring are occupied by Cr'^+ ions, and 
site z = 8 is filled by the M dication. The first term 
is the isotropic, nearest-neighbor Heisenberg-Dirac-Van 
Vleck exchange interaction, with Jcr-Cr and Jcr-M be- 
ing the exchange integrals for Cr-Cr and Cr-M pairs, re- 
spectively. The second term in Eq. ^ describes the local 
crystal-field interaction, whereas the third term gives the 
dipole-dipole intra-cluster interaction, evaluated within 
the point-dipole approximation. Assuming the z-axis 
along the perpendicular to the ring plane, the second- 
order local anisotropy is expected to be dominated by 
the axial terms di, with smaller rhombic terms ef. 
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i=l 



With Si < 3/2, no fourth-order local anisotropy must be 
considered. 

As pointed out by Waldmann [23|, the INS cross- 
section described in for polycrystaUine samples of 
molecular nanomagnets is not of general validity. In par- 
ticular, it may fail to describe correctly the scattering 
from systems with large magnetic anisotropy and does 
not properly take into account intramolecular interfer- 



ence effects. Here we use the following expression, which 
is obtained by averaging the magnetic dipole INS cross- 
section with respect to the possible orientations of 
the scattering vector Q pol |: 
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(3) 

where A — 0.29 barn, Nm is the number of magnetic ions, 
Z is the partition function, ky and ko are the wavevectors 
of the scattered and incident neutrons, exp(— 2M^) is the 
Debye- Waller factor, Q = ko— ky is the scattering vector, 
En is the energy of the generic spin state \n), and the 
function Inm(Q) is defined as: 
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InrniQ) = F* (Q) F, {Q) X \ -MQR^,) + Cln{QR. 
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where F{Q) is the magnetic form factor, Ry gives the 
relative position of ions i and 
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{n\sa,\m){m\s^.\n) {a,j = x,y,z). (6) 



Eq. 0]is the explicit form of the formula given by Wald- 
mann in '2^ . Besides being valid whatever the symmetry 
and amplitude of the magnetic anisotropy, Eq. 01 can be 
easily implemented in a numeric code. 

The parameters of the spin Hamiltonian (Eqs. 1-2) 
have been determined from a best fit of the neutron spec- 
tra based on calculations of Eq. ^ We recall that from 
the INS selection rules only transitions with AS* = 0, ±1 



and AM = 0, ±1 have non-zero intensity. In the fitting 
procedure, a Gaussian lineshape is associated with each 
allowed transition, with a Full Width at Half Maximum 
(FWHM) equal to the instrumental resolution and an 
area proportional to the calculated transition strength. 
For high-energy transitions, a broadening of the final 
state reflecting lifetimes effects has been assumed. The 
parameters defining the Hamiltonian have then been var- 
ied until the best agreement between calculated and ex- 
perimental spectra was obtained. 

The anisotropic part of the spin Hamiltonian can be 
written in terms of irreducible tensor operators (ITO) 
of rank k = 2, and therefore mixes states with differ- 
ent S and M, or at least states with different S if the 
anisotropy is purely axial. Therefore the total spin S is 
not a good quantum number and the total Hamiltonian 
cannot be diagonalized within each (25* -I- l)-dimensional 
block. This difficulty can be overcome by the procedure 
proposed in '2 3] and used to evaluate the mixing between 
the lowest lying S multiplets in the Crg ring |15|. 

First, the minimum ensemble of spin manifolds re- 
quired to reproduce the INS cross-section at high-energy 
transfer is determined assuming isotropic exchange only. 
Then the complete Hamiltonian is diagonalized in the 
corresponding subspace. In this way one obtains the en- 
ergy spectrum and the spin states \n) as linear combina- 
tions of basis vectors \{S)SM) labelled by the set of in- 
termediate spin states (S), with coefficients {{S)SM\n). 

A stringent test for the spectroscopic assignment of 
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the observed transitions is provided by the Q dependence 
of their intensity, which is essentiaUy determined by the 
geometry of the cluster and the composition of the spin 
wavefunctions. This dependence can be easily measured 
with properly calibrated detectors at different scattering 
angles, and the results can be compared with calculations 
using Eq. ^ 



plet with easy-plane zero-field splitting. The INS cross- 
section has been calculated using the expression given 
by Eq. 01 integrated over the Q interval corresponding 
to the experimental conditions. Although quantitatively 
small, S-mixing must be taken into account to reproduce 
correctly the intensity ratio of the observed doublet (the 
intensity calculated with S-mixing neglected is shown in 
Figure 121 by the dashed line). 



IV. EXPERIMENTAL RESULTS 
A. CryNi 

Figure (21 shows the angle-integrated INS intensity 
recorded at T = 2 K for the for CrrNi sample, with IN5 
operated at A = 5 A. The two peaks emerging at 1.19 
meV and 1.34 meV correspond to transitions between the 
\S = 1/2, Af = ±1/2} ground state and the anisotropy 
split sublevels of the \S = 3/2) first excited spin state, 
IS* = 3/2, M = ±1/2} and |S' = 3/2, A/ = ±3/2). 
The solid line represents the intensity corresponding to 
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FIG. 2: Low energy transfer INS response for CrrNi at 2 K (A 
= 5 A). The displayed intensity is the sum for the whole detec- 
tor bank. Background from the sample holder has been sub- 
tracted. Solid line; intensity calculated for the model Hamil- 
tonian described in the text, taking into account S-mixing; 
broken line: calculated response with S-mixing set to zero. 
Inset: intensity calculated with the approximate formula for 
the cross-section given in ref. ll9l| (solid line). 

the Hamiltonian given by Eqs. (1-2), with parameters 
Jcr-Cr = 1-46 meV, Jcr-Ni = 1-69 meV, dcr = —0.03 
meV, and d^vi = —0.6 meV. The dcr value is the one 
determined for the Crg ring and was kept fixed in 
the fitting procedure. As the INS response appears to 
be quite insensitive to the in-plane anisotropy, the non- 
axial part of the Hamiltonian has been neglected. With 
the parameters determined in this work, the ground state 
composition is dominated by the \S = 1/2) component, 
with a small \S = 3/2) contamination (about 1%); the 
first excited state is an almost pure \S = 3/2) multi- 
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FIG. 3: CrrNi INS intensity of the peaks at 1.19 meV (cir- 
cles) and 1.34 meV (triangles) as a function of the scattering 
vector amplitude, Q. Data have been obtained with incident 
wavelength A = 5 A and sample temperature T = 2 K. Cal- 
culated curves are represented by solid lines (to be compared 
with filled circles) and dashed lines (to be compared with 
triangles). Panel a) shows the result obtained with the INS 
cross-section as reported in [l9l |. panel b) shows the output 
of Eq. 01 In both cases Jcr-Cr ~ 1-46 meV, Jcr-Ni ~ 1-69 
meV, dcr ~ —0.03 meV, and djvi — —0.6 meV. The inset in 
panel b) shows the INS response integrated from 1.1 and 1.4 
meV as a function of Q, compared with the calculated curve. 



In the inset of Figure|21 the experimental data are com- 
pared with the curve calculated by the approximate INS 
cross-section reported in ref. (IQj . and the Hamiltonian 
parameters quoted above. Due to the large anisotropy at 
the Ni site, the Borras-Almenar et al. approximation is 
not adequate in the present case, and its use instead of 
Eq. 01 would lead to a completely wrong estimate of the 
Hamiltonian parameters. This difference is more evident 
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when comparing the Q dependence of the transition in- 
tensities with theoretical estimates provided by the two 
formulae. As shown in Figure IStj the intensity of the 
transition involving the \S = 3/2, M — ±1/2) excited 
level at 1.19 meV is underestimated over a large Q inter- 
val by the expression given in |19| | , whilst the intensity of 
the transition towards the \S = 3/2, M = ±3/2) state at 
1.34 meV is overestimated by more than a factor two. In 
comparison, an excellent agreement between calculated 
curves and observed data is obtained when using Eq. ^ 
(Figure OJj) . 




Energy Transfer (meV) 



FIG. 4: CrrNi INS spectra collected with A = 4A atT = 2 
K (circles) and T = 12 K (triangles). Smooth lines represent 
the spectra calculated from eigenvalues and eigenvectors of 
the Hamiltonian given in Eas lll2l by associating to each tran- 
sition a Gaussian lineshape with a height proportional to the 
calculated probability and a width of 147 peV, corresponding 
to the experimental resolution. 



Excitations involving spin manifolds with higher en- 
ergy can be observed by reducing the incident wave- 
length. Results obtained with A = 4A atT=:2K 
and T = 12 K are shown in Figure H At T = 2 K, 
weak excitations are observed between 3 and 4 meV. 
They correspond to transitions from the ground state 
to \S = 1/2) and 15* = 3/2) excited states. When the 
sample is warmed at 12 K, a strong peak appears at 2.08 
meV. This excitation is attributed to the transition from 
the first excited \S = 3/2) state to a jS* = 5/2) state lying 
at 3.31 meV. The calculated spectra, represented by the 
smooth lines in Figure ^ compare very favorably with 
the experiment and confirm the validity of the proposed 
set of parameters. The ratio Jcr-Ni/ Jcr-Cr = 1-16 is 
somewhat larger than the recent estimate based on the 
fit of heat capacity and torque magnetometry data |0| . 
In addition the weighted average axial anisotropy param- 
eter d = —0.045 meV 0| is about 1.7 times larger than 
previously reported, because of the very large anisotropy 
of the Ni ion that was not estimated in previous works. 



B. CrrZn 

Substitution of a diamagnetic Zn^+ ion for one Cr'^^ 
in the AF Crg ring results in a jS* = 3/2) spin ground 
state, whose anisotropy splitting can be observed by high- 
resolution INS experiments. Figure [S] shows the spectra 
recorded on IN5 with incident wavelength 9 A (energy 
resolution at the elastic peak of 19 /ieV), with the sample 
kept at 2 and 10 K. 

3 I r- ^ ■ ^ ■ 1 
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FIG. 5: CrrZn high-resolution INS spectra collected with A 
= 9A atr = 2K (circles) and T = 10 K (triangles). The 
peak at 0.11 meV is an intra- multiplet transition involving the 
anisotropy split components of the |S = 3/2 > ground state. 
Solid (2 K) and dashed (10 K) lines are intensities calculated 
for a purely axial, easy plane magnetic anisotropy of the Cr 
ions {dcr = —0.028 meV). The elastic peak and a quasielastic 
contribution have been included in the best fit procedure. 

Assuming crystal-field parameters similar to those de- 
termined in the Crg parent compound, we expect an 
easy- axis anisotropy for the S — 3/2 ground manifold 
and an easy-plane splitting for the lowest 5 = 5/2 multi- 
plet. In this hypothesis, the peak observed at about 0.11 
meV must be attributed to the \S = 3/2, M = ±3/2 >-> 
\S — 3/2, M = ±1/2 > transition, as only the lowest 
manifold is thermally populated at T = 2 K. The solid 
and broken lines give the calculated intensity for T — 2 
and 10 K assuming a purely axial local crystal field for 
the Cr ions, with dcr = —0.028 meV [dzn =0). As for 
the Ni case, the INS spectra are not sensitive to the in- 
plane component of the anisotropy and again we consider 
only the axial part of H . 

As shown in Figure EJ intermultiplet transitions at 
higher energy transfer are observed with incident wave- 
lengths of 4 and 5 A. The peaks appearing in the T = 2 K 
spectrum are due to transitions from the split \S = 3/2 > 
ground state towards \S = 1/2 > (0.83 meV), \S = 5/2 > 
(1.91 meV), IS" = 1/2 > and \S = 3/2 > (2.25 meV), 
IS" = 1/2 > (3.0 meV), and unresolved jS" = 5/2 > and 
\S = 3/2 > states (3.6 meV). 

The spurious peak at 0.4 meV in panel b is due to neu- 
trons diffused incoherently by the sample and diffracted 
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FIG. 6: INS spectra for CryZn at 2 K (circles) and 12 K (tri- 
angles), measured with incident wavelength of (a) 5 A and 
(b) 4 A . The 0.4 meV peak in panel b is spurious. Smooth 
lines are best fits of the experimental data to a superposi- 
tion of Gaussian line-shapes associated to allowed transitions. 
Each Gaussian contribution has an area proportional to the 
transition strength calculated with Eq. 0] assuming isotropic 
exchange {Jcr-Cr = 1.43 meV, Jcr-Zn = 0) and axial local 
anisotropy {dcr ~ —0.028 meV). 



by the cryostat walls. The instrumental resolution in this 
configurations is not high enough to resolve the small 
anisotropy splitting of the excited levels, so that the en- 
ergy separation of the doublets around 0.83 and 1.91 meV 
reflects the splitting of the ground state. Hot peaks as- 
sociated to excitations from the first excited jS* = 1/2 > 
level appear at T = 12 K around 1.39 and 2.6 meV. 
The solid and broken lines in Figure are intensities cal- 
culated at the two temperatures assuming an exchange 
parameter v ery similar to the one reported for the Crg 
ring m ref. [la, Jcr-Cr = 1.43 meV for every nearest- 
neighbor Cr-Cr pairs, Jcr-Zn — (as Zn is diamagnetic) 
and dcr — ~0.028 meV. The agreement with experimen- 
tal data is very good and justifies the use of just one value 
for the Jcr-Cr exchange integrals, although the use of a 
less symmetric magnetic model would improve the fitting 
at high energy transfer. The spectroscopic assignments 
are confirmed by the Q dependence of the transition in- 
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Momentum Transfer (A"') 

FIG. 7: Q dependence of the peaks appearing in the CrrZn 
inelastic spectra shown in Figure|^ compared with theoretical 
estimates (solid lines). Triangles \S = 3/2 >^ \S = 1/2 > 
at 0.83 meV; circles 15 = 1/2 >^ \S = 3/2 > at 1.38 meV; 
squares IS' = 3/2 >— > \S = 5/2 > at 1.91 meV; open circles 
\S = 3/2 >^ \S = 5/2 > at 3.6 meV. 



tensities, as shown in Figure |3 where experimental data 
for some of the excitations observed are compared with 
curves calculated with Eq. 01 It is interesting to notice 
that the general behavior predicted for cyclic spin clus- 
ters by Waldmann 20] is followed also in the present case 
where the ring symmetry is broken by the dication. 



C. CrrMn 

Substitution of one high-spin, s = 5/2, Mn^+ ion for 
a Cr'^"'" in the Crg ring results in an uncompensated 
IS* = 1 > ground state. In the presence of purely axial 
anisotropy, the zero-field splitting would lead to a quasi- 
triplet structure, with an \S = 1,M — ±1 > doublet 
and an \S = 1,M = > singlet separated by an en- 
ergy proportional to the axial crystal field parameter. As 
shown in Figure|Hl the A = 9 A high-resolution spectrum 
recorded at T = 2 K shows a broad excitation that can be 
fitted by the superposition of two Gaussians, each with 
a FWHM equal to the instrumental resolution (19/ieV) 
and centered at 0.095 and 0.11 meV, respectively. This is 
a clear indication of a sizeable rhombic term in the zero- 
field splitting Hamiltonian. By fixing the dcr parameters 
to the values determined by INS for the parent Crg com- 
pound {dcr = —0.03 meV) 15], the best fit of the data 
is obtained for dMn = —3 fieV (a factor two hundred 
smaller than the Ni case) and eMn — scr — 0.14djv/T!- 
If the S mixing is neglected, the ZFS parameters de- 
scribing the anisotropic splitting of the S* = 1 ground 
multiplet can be obtained by projecting the complete 
Hamiltonian onto the corresponding subspace. We ob- 
tain Hs^i = D{Sl - S{S + l)/3) + E{Sl - SD with 
a rhombicity E/D ~ 0.13. The latter is reduced by S 
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FIG. 8: CrrMn high-resolution INS spectra collected with A 
= 9 A at T = 2 K. Two intra-multiplet transitions within 
the anisotropy split components of the |5 = 1 > ground 
state are observed. The solid line is the intensity calcu- 
lated assuming both axial and in-plane anisotropy. The dash- 
dot lines represent the individual contributions to the broad, 
unresolved excitation from the transitions between the split 
IS = 1,M = ±1 > quasi-doublet and the |S = 1,M = > 
singlet. The dot line gives the intensity calculated with S- 
mixing set to zero. The elastic peak and a quasielastic con- 
tribution have been included. 



mixing to E/D = 0.1 (see Fig. IS)). 

Intermultiplet excitations measured with incident 
wavelength A = 4 A are shown in Figure El At T = 2 K 
only one strong magnetic peak is observed at 1.6 meV. 
A second magnetic transition appears at 2.3 meV if the 
sample is heated at 12 K. The two peaks correspond to 
the transitions involving the ground state jS* = 1 > and 
the \S = 2 > and jS* = 3 > multiplets with minimal 
energy. 

Both peaks therefore correspond to excitations within 
the lowest rotational-like band and their energy follows 
the Lande rule very closely. In these conditions, only an 
average exchange parameter can be estimated. Assuming 
Jcr-Cr = Jcr-Mn = 1.43 mcV, Calculations provide the 
spectra shown in Figure El by smooth lines. The agree- 
ment is satisfactory and a model with more parameters 
would not be justified by the data. Alternatively, by fix- 
ing Jcr-Cr — 1.46 meV as in CryNi, the same results 
are obtained with Jcr-Mn = 1.37 meV. The assignment 
of the two transitions to the Lande-band is corroborated 
by the Q-dependence of their intensity, shown in Figure 
[Tni For both peaks a pronounced oscillatory behavior is 
observed, whilst an almost flat Q dependence is expected 
for transitions to states not belonging to the rotational 
band. 
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FIG. 9: Inelastic spectra at T = 2 K (circles) and T = 12 K 
(triangles) obtained for CryMn with an incident wavelength 
A = 4 A. The peak at 0.4 meV is due to spurious effects. 
Smooth lines are intensities calculated at 2 K (solid) and 12 
K (dashed). The elastic peak and a quasielastic contribution 
have been included. 
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FIG. 10: Q-dependence of the peaks appearing in the CrTMn 
inelastic spectrum at 12 K, compared with theoretical esti- 
mates: |S'=1>— >|5' = 2>atl.6 meV (triangles and solid 
line); |S = 2 >^ 15* = 3 > at 2.3 meV (circles and dashed 
line). The Q-dependence of the intramultiplet transition at 
0.1 meV (T = 2 K) is compared with calculations in the inset. 



V. DISCUSSION 

The energy spectra of the heterometallic Cr7M mag- 
netic wheels resulting from our INS investigation are 
shown in Figure [TTl as functions of the total spin, S. A 
parabolic band, formed by states with minimal energy 
for each S value, is easily identified. The levels belonging 
to this parabolic band have energies that closely follow 
the Lande interval rule, Es = Aio[5(5 + 1) - So{So + 
l)]/[5i(S'i + 1) - So{So + 1)], where 5*0 is the spin of 
the ground state and Aio is the energy of the first ex- 
cited level, with spin ^i. Excitations involving adjacent 
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a crossing with the ground state at B = 6.9 T with no 
level repulsion and a vanishingly small Schottky anomaly 
in the heat capacity On the other hand, according 
to the results of the diagonalization of the Hamiltonian 
Eqs. I1I2I anticrossings between the lowest lying levels are 
expected in the CrrM compounds. At the anticrossing 
conditions, a finite Schottky anomaly will occour, with 
an amplitude that depends on the angle between the ap- 
plied magnetic field and the easy-axis. Moreover, the 
enhancement of S'-mixing at the level anticrossings will 
produce maxima in the torque signal corresponding to 
oscillations of the total spin amplitude, as observed for 
the Mn[3 x 3] grid [113. 

In the investigated series, CryNi is particularly inter- 
esting. The amount of S'-mixing in this ring is quite 
small, corresponding to about 1% admixture of S ^ 1/2 
components in the ground state. CryNi can therefore be 
considered as an effective 5 = 1/2 system suitable for the 
implementation of the qubit [13|| . 



FIG. 11: Energy of the lowest spin eigenstates as a function 
of the total spin S calculated for CryM magnetic wheels; a) 
M = Cr, b) M = Ni, c) M = Mn, d) M = Zn. Isotropic 
exchange interactions only are assumed, with parameters a) 
Jcr-Cr = 1.46 meV b) Jcr-Cr =1.46 meV and Jcr-Ni 
= 1.69 meV, c) Jcr-Cr =1.43 meV and Jcr-Mn ~ 1-43 meV, 
and d) Jcr-Cr =1.43 meV and Jcr-Zn = 0. 



levels of the parabolic band have intensities with similar 
Q dependencies, with an oscillatory behavior and a pro- 
nounced maximum at a Q value related to the radius of 
the wheel Q . As discussed in 0, , these excitations are 
related to the combined quantum rotation of the oppo- 
sitely oriented total spin on each Neel sublattice of the 
AF wheel |25|. 

Effects due to the mixing of different spin multiplets 
have been considered. Such effects proved to be impor- 
tant to correctly reproduce the energy and intensity of 
magnetic excitations in the neutron spectra. An inter- 
esting difference between the spin wavefunctions of the 
parent Crg ring and those of the heterometallic deriva- 
tives concerns their symmetry. In Crg the ground state 
and the first excited state belong to different irreducible 
representations of the molecule point group, whereas they 
have the same symmetry in the substituted wheels. This 
has important consequences on the system behavior in 
the presence of an external magnetic field. In Crg the 
Zeeman splitting of the \S = 1,±1 > doublet leads to 



VI. SUMMARY 

Intra-multiplet and inter-multiplet excitations involv- 
ing spin manifolds with energy smaller than 4 meV have 
been measured with Inelastic Neutron Scattering in poly- 
crystalline samples of the heterometallic wheels CrrM 
(M = Ni, Mn, Zn). The minimum set of exchange and 
local crystal field parameters necessary to describe the 
physics of each investigated compound has been deter- 
mined by comparing the experimental spectra with the- 
oretical cross-sections. The results obtained show that 
chemical substitution of magnetic ions in a cyclic struc- 
ture can be used to tailor the magnetic properties of the 
wheels, by controlling the microscopic exchange inter- 
action. Finally, the explicit form of the general powder- 
averaged INS cross-section given in Ref. j23 and suitable 
to describe molecular nanomagnets of any symmetry has 
been presented. 
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